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Abstract 

We explore numerically the feasibility of using chemical patterning to control the size and polydispersity 
of micron-scale drops. The simulations suggest that it is possible to sort drops by size or wetting properties 
by using an array of hydrophilic stripes of different widths. We also demonstrate that monodisperse drops 
can be generated by exploiting the pinning of a drop on a hydrophilic stripe. Our results follow from using 
a lattice Boltzmann algorithm to solve the hydrodynamic equations of motion of the drops and demonstrate 
the applicability of this approach as a design tool for micofluidic devices with chemically patterned surfaces. 
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I. INTRODUCTION 



In recent years advances in patterning techniques, as well as in our understanding of electrowet- 
ting, have made it increasingly feasible to construct surfaces with regions of different wettability 
on micron length scales. The behaviour of fluids spreading and moving across such surfaces is 



extremely rich, and is only just beginning to be explored ll 



Mm 




Id, 



ilD. Biosys- 



tems have evolved to use hydrophobic and hydrophilic patches to direct the motion of fluids at 



surfaces [tl2|, 



130 . Examples are desert beetles, whose patterned backs help them to collect dew 



from the wind, and the leaves of many plants, designed to aid run-off of rain water. Similar surface 
patterning have also been exploited in the design of microfluidic devices. For example, Zhao et al 



lll4 [150 have used chemically patterned stripes inside microchannels to control the behaviour of 
liquid streams and Handique et al [fl6] have used a hydrophobic surface treatment to construct 
a device which can deliver a fixed volume of fluid. Our aim here is to use numerical simulations 
to further demonstrate ways in which chemically patterned surfaces can be used to manipulate 
the behaviour of micron-scale drops, demonstrating that they provide a useful tool for exploring 
device designs. 

One application where one might envisage a use for surface patterning is in controlling the size 
and polydisper sity of an array of drops, an important consideration in many microfluidic technolo- 



gies, e.g. mmV 
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and the references therein. We 



explore numerically two ways in which this might be possible. The first uses hydrophilic stripes 
of different widths to sort drops by size or wetting properties. The second produces monodisperse 
drops by exploiting their tendency to remain pinned on hydrophilic regions of the surface. 

Most drop sorting techniques available in the literature involve active actuation to control the 
motion of the drops. This can be done by (but is not limited to) electro-osmotic [21], mechanical 
P22I] . optical [123D, or dielectrophoretic manipulation. Our proposed devices are passive de- 
vices, which require neither a detection nor switching mechanism. Tan et al. [25] make use of 
the channel geometry to control the flow and hence sort drops by their size. Here we propose a 
different drop sorting device which uses chemical surface patterning to separate the liquid drops. 
Since a given surface is wet differently by different liquids, our proposed method can be used to 
separate drops based on their wetti ng p roperties, size, or both. 

271] and passive [128|, HQ] drop generating techniques are known in the literature. 



Active [|26, 



but to the best of our knowledge, they are mostly concerned with generating liquid drops from 
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continuous streams or larger sample drops. Here, we look at the problem from a different angle 
and suggest how monodisperse liquid drops can be formed from small poly disperse drops. 

Since capillary phenomena are increasingly important for systems with large surface to volume 
ratio, passive microfluidic devices based on surface chemical patterning can play an important role 
in the future, particularly as the devices are downscaled. Careful designs of capillary networks may 
even lead to an integrated device, where drop generation, sorting, mixing, and delivery can all be 
done on a single chip. We further note that it is an intriguing idea to incorporate electrowetting into 
both devices, adding active control, flexibilities, and functionalities to the devices. Such devices 
are not limited to the liquid- gas system considered in this paper. The designs can be used for both 
a liquid- liquid and liquid-gas system, and in both an open and closed channel. 



We use a lattice Boltzmann algorithm [|31 



330 to solve the hydrodynamic equations of 



motion for the drop movement across the patterned substrates thus exploring the key parameters 
of possible devices. This numerical approach has been shown to agree well with experiments 
in several cases 111, l2l, giving us confidence in using it for the more complicated situations 
considered here. 



II. EQUATIONS OF MOTIONS 

The equilibrium properties of the drop are described by a continuum free energy 

^ = f {Mn) + ^{danY)dV+ I ^s{ns)dS. (1) 
Jv 2 Js 

ipb{^) is a bulk free energy term which we take to be [13 if] 

Mn) = Pc{^n + lf{^l - 2zy, + 3 - 2/3r^) , (2) 

where z/^ = (n — n^juc, = (T^ — T)/Tc and n, ric, T, and Pc are the local density, crit- 
ical density, local temperature, critical temperature and critical pressure of the fluid respectively. 
/? is a constant typically chosen to be 0.1. This choice of free energy leads to two coexisting 
bulk phases of density nc(l ± ^/W^)- The second term in Eq. ([T]) models the free energy as- 
sociated with any interfaces in the system. is related to the liquid-gas surface tension via 
(Jig = (4y2^p^(/3r^)^/^nc)/3 till]. ^Th e last term describes the interactions between the fluid and 
the solid surface. Following Cahn [34] the surface energy density is taken to be V^s(n) = —(pUg, 
where Ug is the value of the fluid density at the surface, so that the strength of interaction, and 
hence the contact angle, is parameterised by the variable 0. 
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The dynamics of the drop is described by the continuity and the Navier-Stokes equations 

dtu + da{nua) = , (3) 

dt{nUa) + df3{nUaUf3) -df3Pa(3 + na^ + iydf3[n{df3Ua + daUf3 + 5af3djU^)] , (4) 

where u, P, i/, and a are the local velocity, pressure tensor, kinematic viscosity, and acceleration 
respectively. We impose no-slip boundary conditions u = on the surfaces. 

The thermodynamic properties of the drop are input via the pressure tensor P which can be 
calculated from the free energy 

Paf3 = {Ph - ^{do^nf - Knd^^n)5af3 + i^{do,n){d^n) , (5) 
where = Pc{^n + ^f[^^l - 2z/^ + 1 - 2/?r^) and /x^ = ^n^l^h\n=nk = - ^^^). 



III. LATTICE BOLTZMANN 



The lattice Boltzmann algorithm is defined in terms of the dynamics of a set of real numbers 
which move on a lattice in a discrete time. A set of distribution functions {/i(r, t)}, is defined 
on each lattice site r. Each of these distribution functions can be interpreted as the density of 
the fluid at time t that will move in direction i. The directions {i} are discrete and for a three 
dimensional system, one needs to take at least fifteen velocity vectors = c(0, 0, 0), c(±l, 0, 0), 
c(0, ±1, 0), c(0, 0, ±1), c(±l, ±1, ±1). c = Ax/ At is the lattice speed, and Ax and At represent 
the discretisation in space and time respectively. The distribution functions are related to the 
physical variables, the fluid density n and its momentum nu, through 

Y^fi = ^, fi^ia = nUa. (6) 

i i 

Taking a single-time relaxation approximation, the evolution equation for a given distribution 
function fi takes the form 

Mr + ViAt, t + At) = Mr, t) + /''(^^ ^) " ^) + nw^v^^g^ . (7) 

r 

where = wi = 1/3 if |^| c and = W2 = 1/24 if l^'l = c\/3. The relaxation time r 



tunes the kinematic viscosity via u = ((Ax)^(r — l/2))/(3At) 113 111 . The last term in Eq. (|7]) 



is the forcing term, which is related to the acceleration through = Since Ax and At 
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are typically taken to be 1 in simulation units, the notations and ga are interchangeable in this 
paper. 

It can be shown that Eq. (|7]) reproduces Eqs. dS]) and d?]) in the continuum limit if the correct 
thermodynamic and hydrodynamic information is input to the simulation by a suitable choice of 
local equilibrium functions, i.e. if the following constraints are satisfied 

E /r = n , E /r^ia = nu^ , (8) 

i i 
i 

r.2 



i ^ 

A possible choice for f^^ is a power series expansion in the velocity 
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G2aa = 



Full details of the lattice Boltzmann algorithm are given in 



infl 



33f] 



(9) 



IV. SORTING DROPS BY SIZE AND WETTING PROPERTIES 

To use chemical patterning to sort drops according to size or wettability we consider the design 
in Fig. 1. A surface is patterned with a rectangular grid of hy drophilic (relative to the background) 
stripes. A drop is input to the device at A and subject to a body force at an angle < 45"" to the 
X-axis. The system is confined in a channel of height L^. 

The path taken by the drop through the device depends on the drop contact angles with the 
substrate and the strength of the body force. It also, of particular relevance to us here, depends on 
the width of the stripes relative to the drop radius. By choosing the stripes along the y direction 
to be of equal widths, but those along x to increase in width with increasing drops of different 
sizes move along different paths with the larger drops moving further along the y direction. 

Fig. 2(a) - (c) show simulations of the paths of drops of initial radius R = 25, 26 and 29 
moving through such a device. The simulation parameters are: 6^phiiic = 60^, 6^phobic = 110^, 
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= 80, = 3.0 10-^ ay = 2.0 10"^ dynamic viscosity r] = 0.4128, aig = 7.710-^ drop 
density ni = 4.128, Si = 20, S2 = 30, S3 = 40, Sy = 20, Sh = 80, and Ss = 200. The drop 
velocity is measured to be ^ 5 . 10~^. Comparing the simulation results with experiments for 
a simpler surface pattern [jl], we established a mapping between simulation and physical units, 
except that the drops move too fast by a factor of ^ 500 in the simulations. (This is a consequence 
of the interface width being too wide relative to the drop size in this, and indeed all, mesoscale 
simulations.) Taking into account this empirical factor, the parameters used here correspond to 
dimensionless numbers: Re = riivR/r] ^ 3 10~^, Ca = rjv/aig ^ 6 10~^, Bo = niaR^ /aig ^ 
1.2, and (3^ = niv^R/aig - 2 10"^ 

When a submilimetric drop of initial radius R is jetted on a flat homogeneous substrate, it will 
relax to form a spherical cap with a contact angle given by the Young's law. This is however not 
generally the case on chemically patterned surfaces. The equilibrium drop shape at A is elongated 
in the x-direction because the drop prefers to wet the hydrophilic stripe. The introduction of the 
body force will further deform the drop shape so that it is no longer symmetric in either the x 
or the ^-directions. Indeed, the drop will be confined in a hydrophilic stripe only if the effective 
capillary force is able to counterbalance the imposed body force. 

In cases where the drops are confined in the Si stripe, they will move in the x-direction from 
A to the cross-junction where their paths may diverge. In order for a drop to move in the y- 
direction, the capillary force in this direction must be large enough to overcome the sum of the 
capillary force and the excess external body force in the x-direction (recall ax > ay). This is where 
the asymmetry of the drop shape comes into play. As the volume of the drop is increased, a larger 
fraction of it overhangs the stripes and hence a larger fraction will interact with the hydrophilic 
stripe along the y-direction at the junction. This increases the capillary force along y and means 
that larger drops (e.g. R = 26) will move in the ^-direction to point C, whereas smaller drops (e.g. 
R = 25) will continue to move along x. 

Since ax > ay and S2 > Sy the asymmetry of the drops' shape guarantees that drops at C move 
to D. If the widths of the stripes along x were all the same the drops would then move from D to 
E for the same reason that they moved from B to C. However, increasing the stripe width so that 
S2 > Si reduces the effective capillary force in the y-direction and therefore some of the drops 
(e.g. R = 26) continue to move along x whereas larger drops are routed around the corner towards 
the third stripe. Hence the drops are sorted by size. 

Fig. 2(c) shows the path of a drop of radius R = 29 through the device. Although it does 
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move up to the third stripe it has to pass several junctions before it does so. This occurs because 
the drop takes time to relax to its steady state shape as it moves along the 82 stripe - note that 
the drop morphology is slightly different as it crosses the second and third vertical stripes, with 
a slightly increased overhang at the latter. Ideally the distance between two vertical hydrophilic 
stripes should have been longer to overcome this effect. 

These simulations suggest that by increasing the number of stripes and carefully controlling 
their widths it may be possible to sort polydisperse drops into collections of monodisperse drops. 
The precision with which this can be done is determined by the increase in the stripe widths. Here 
the vertical stripe width is fixed. As a result {Rn+i — Rn) > {Rn — Rn-i)- {Rn+i — Rn) ^ 
{Rn — Rji-i) can be achieved by increasing the vertical stripe width with increasing x. Two other 
parameters, the wettability contrast and the external body force, could also be adjusted to fine-tune 
the device. 

This design can also be applied to sort drops of similar size that possess different wetting 
properties. Returning to Fig. 1 all the drops will move to the right from point A to point B. Drops 
with a lower contact angle on the stripes will overhang the stripes less. Hence at B, these drops will 
interact less with the stripes along y and will continue to move along x. As the contact angle of 
the drops increases they will feel more effect from the junction and, at some threshold wettability, 
they will turn the comer to move along y. As shown in Fig. 3, we find that for drops of radius 
i? = 26 the drop moves along stripe 5i for ^phiuc = 40° but is diverted to stripe 82 for 6^phiiic = 60°, 
with other simulation parameters as before. 

The current device speed and throughput are limited by the need to eliminate unwanted drop 
coalescence. For example, in the drop sorter, drops of different sizes and wetting properties move 
at a different speed. If two drops are entered too closely together, it may lead to a situation where 
the two drops coalesce and mix instead of being sorted. 

V. GENERATING MONODISPERSE DROPS 

As a second example we demonstrate how suitable chemical patterning can be used to generate 
monodisperse drops. In Fig. 4(a) regions A and C are hydrophobic to the drop relative to regions B 
and D. For simplicity, we take = = 110° and Ob = 0^ = 60°. Polydisperse drops are input 
to the system in region A and subject to an external body force in the x-direction. If the applied 
body force is large, all drops will exit the system at region D without change. 
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If the body force is small, however, the effective capillary force at the B-C border can balance 
the external force and the drops are trapped in the hydrophilic stripe. As more and more poly- 
disperse drops coalesce at the B-C border, the resultant drop will no longer be confined to the 
hydrophilic stripe. A fraction of the drop volume will start to occupy the hydrophobic region. At 
a certain critical volume, the front end of the drop will reach the D region and the drop will tunnel 
from B to D. This is shown in Fig. 4 where we have presented the lattice Boltzmann simulation 
results for the following set of parameters: 5a = 60, 5b = 70, 5c = 30, 5d = 140, and = 10~^. 
Other parameters remain the same. 

The accuracy with which it is possible to control the final drop size is dependent on the ratio 
between the average input drop volume to that of the generated output drop: the smaller the ratio 
the better the accuracy. There are three key parameters that can be used to tune the generated drop 
volume: the wettability contrast, the stripe widths, and the external body force. Increasing either 
the wettability contrast or the stripe widths will increase the effective capillary force, and hence 
increase the final drop volume. Increasing the body force will have the opposite effect, resulting in 
a smaller output drop. It is also important to tune the input drop frequency. If the frequency is too 
high then a number of input drops will coalesce with the resultant drop as it moves from region B 
to D. 

VL DISCUSSION 

To conclude, we have proposed two ways in which chemically patterned surfaces can be used 
to control drop size and polydispersity. The drop sorter takes a collection of polydiperse drops as 
input and separates them into different channels based on their size or wetting properties. The drop 
generator allows small polydisperse drops to coalesce, pinning the resultant drop until it reaches a 
critical volume. 

The behaviour of drops within the devices was studied by solving their equations of motion 
using a lattice Boltzmann algorithm. Thermodynamic properties of the drops, such as surface 
tension and contact angles were included by using a simple free energy model: in equilibrium the 
drops minimise a Landau free energy functional. 

Our results show the potential of chemical patterning as a way to control the behaviour of 
micron-scale drops and demonstrate that simulations provide a useful tool for exploring possible 
device designs. Further progress needs experimental input to help to understand more fully the 
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mapping between physical and simulation parameters and the experimental constraints in the re- 
alisation of these devices. We note that the simulations are computer-intensive, for example the 
results in Fig. 2(c) took about a month on 8 nodes of dual 2.8 GHz Xeon machines. Synergy 
with experiment will help focus computer resources on the most realistic and useful questions and 
parameter ranges. 
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FIGURE CAPTIONS 

Fig. 1: Schematic diagram of a drop sorter. The grey stripes on the surface are hydrophiUc with 
respect to the background. 5 labels the widths of the stripes and a the imposed acceleration. The 
arrows show possible paths of a drop through the device. 

Fig. 2: Paths taken by drops of radius (a) R = 25, (b) R = 26, and (c) i? = 29 through the drop 
sorter. 5i = 20, 82 = 30, ^3 = 40, and Sy = 20. 

Fig. 3: Paths taken by drops of radius i? = 26 through the drop sorter, when the equilibrium 
contact angle of the drop on the hydrophilic stripes is (a) 6^phiiic = 40° and (b) ^phiiic = 60°. 

Fig. 4: Simulations showing the evolution of the drop shape and position with time in the drop 
generator, (a) Contour plot as volume is added to the system, (b)-(e) Snapshots as volume is 
added, (b) V = 1.2 10^ (c) V = 2.6 10^ (d) V = 3.3 10^ (e) V = 4.1 10^ 
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FIG. 3: 
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